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Abstract 

Ml transitions from the ^Li(0+;T = 1) state at 3.563 MeV to the •^Li(l+) ground 
state and to the a + d continuum are studied in a three-body model. The bound 
states are described as an a + n + p system in hyper spherical coordinates on a 
Lagrange mesh. The ground-state magnetic moment and the gamma width of the 
^Li(0"^) resonance are well reproduced. The halo-like structure of the ^Li(0"^) res- 
onance is confirmed and is probed by the Ml transition probability to the a + d 
continuum. The spectrum is sensitive to the description of the a+d phase shifts. The 
corresponding gamma width is around 1.0 meV, with optimal potentials. Charge 
symmetry is analyzed through a comparison with the /3-delayed deuteron spec- 
trum of ^He. In ^He, a nearly perfect cancellation effect between short-range and 
halo contributions was found. A similar analysis for the ^Li(0^;r = 1) 7 decay is 
performed; it shows that charge-symmetry breaking at large distances, due to the 
different binding energies and to different charges, reduces this effect. The present 
branching ratio r^(0"'" — > a + d)/r^(0"*" 1^) 1.3 x 10~^ should be observable 
with current experimental facilities. 
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1 Introduction 



Electromagnetic transition processes provide a useful tool for the study of 
the nuclear structure and of the reaction mechanisms. The theoretical study 
of such processes yields estimates for the different static and dynamical ob- 
servables of a nucleus. In ^Li, the {0'^;T = 1) state has raised interest as a 
good candidate for observing parity violation [112] ■ Indeed, its decay into the 
a + d continuum is forbidden by parity conservation. Since electromagnetic 
Ml transitions into this continuum are allowed, they have been also stud- 
ied because they may compete with the parity-violating decay and make its 
detection difficult. 

However, the ^Li(0"'"; T = 1) state is also interesting by itself. It is most likely 
a halo state, as it is the isobaric analog of the ^He ground state [3]. Ml 
transitions to the continuum are in fact also an excellent tool to explore these 
halo properties and compare them with those of ^He. Recent experimental 
and theoretical works on the delayed ^He (3 decay suggest that the deuteron 
spectrum is strongly sensitive to the halo structure (see Ref. |1] and references 
therein). Similarities between this process and the 7-delayed deuteron emission 
of ^Li(O^) are expected, and should test charge symmetry in exotic light nuclei. 

The branching ratio of the total transition probability to the a + d continuum 
and the transition probability to the ^Li(l"'") ground state was estimated as 
8 X 10~^ under a number of simplifying assumptions However the shape 
and magnitude of the transition probability to the continuum as a function 
of the deuteron energy were not studied. In addition, the sensitivity with 
respect to the a + d potential, as well as convergence problems, were not ad- 
dressed. The aim of the present work is to investigate Ml transitions from the 
^Li(0"'") excited state to the a + d continuum, as well as to the ^Li(l"'") ground 
state. For the description of the ^Li states, we use different two-body and 
three-body models. Three-body hyperspherical wave functions [5j are based 
on the Lagrange-mesh method and give an accurate solution of the three-body 
Schrodinger equation [B]- The a + d scattering wave function is factorized into 
a deuteron wave function and a nucleus-nucleus scattering state. This work 
extends our previous study on ^He /5 decay where the same formalism was 
used. 

The model is presented in Section 2. The potentials and the corresponding 
two-body and three-body wave functions are also described. In Section 3, we 
discuss the results in comparison with the experimental data, and analyze the 
sensitivity with respect to the a + d potential. Finally, conclusions are given 
in Section 4. 
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2 Model 



2.1 Three-body wave functions of^Li bound states 



The ^Li bound-state wave functions are defined in an a + n+p model using the 
hyperspherical coordinates [5]. A set of Jacobi coordinates for three particles 
with mass numbers Ai = 1, A2 = I, and ^3 = 4 is defined as 



x = ^/JI^r, y = ./JI^R, (1) 

where the (dimensionless) reduced masses are given by /ip„ = 1/2 and /i^^ = 
4/3. The p — n relative coordinate and the coordinate between a and d are 
denoted by r and R, respectively. Equations ([1]) define six coordinates which 
are transformed to the hyperspherical coordinates as 

= + y'^, a = arctan(?//x), (2) 

where a varies between and it/2. With the angular variables = {dxiVx) 
and VLy = {9y,ipy), equations define a set of hyperspherical coordinates 
which are known to be well adapted to the three-body Schrodinger equation. 

We define 7 = {i^, ^y, L, S) where and iy are the orbital momenta associated 
with the Jacobi coordinates x and y, respectively. With the notation = 
(a, fix, fly), a three-body wave function with spin J and parity tt reads [6] 

^il(''{p,n,) = p~'/'Y.xiKip)y'Km, (3) 

■yK 

where y:^^^{fl5) are the hyperspherical functions (including spin), defined as 



JM 



(j)^^^y (a) = Uk^" (cos aY^ (sin aY^Pn 



(cos2q;), 



(4) 



with n = [K — i-x ~ ^j/)/2, and where J\f^^^ is a normalization factor, 
a Jacobi polynomial and ^ spin function. Details are given in Refs. 



The hyperradial functions x^k cire obtained from a set of coupled equations 
truncated at K = K^ax- The problem is solved by using the Lagrange-mesh 
technique (see Ref. [6] for details). 



The three-body wave functions contain components with total intrinsic spin 
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S = and 5* = 1. Because of the positive parity, ix + iy is even and only even 
K values are involved. 



2.2 a + d two-body wave functions 

As it was done in Ref. [Ij for the ^He /3 decay, the scattering a + d wave 
functions are factorized into a deuteron ground-state wave function, calculated 
with an appropriate NN potential, and an a + d wave function derived from 
a potential model. We neglect the small D component of the deuteron. In the 
a + d exit channel, only S waves are involved. Consequently, the final 1"*" wave 
function reads 

r, R) = vl/,(r) ^UE; R). (5) 
The spatial part of the deuteron wave function is written as 

^d{r) = r-'ud{r)Yoo{r). (6) 



The S'-wave component of the a + d relative motion wave function is factorized 

as 

^adiE;R) = R-^UE{R)YooiR). (7) 



The normalization of the scattering wave function is fixed by the asymptotic 
behaviour as 



ue{R)-^ FoiKdR) cos5o(^) + GoiKdR) sin5o(^), 

it— >oo 



where Fq and Go are the Coulomb functions, So{E) is the s-wave phase shift 
at energy E, and kad is the wave number of the relative motion. 

The present two-body model can also be applied to a + d bound states. In that 
case, the scattering wave function ue{R) in Eq. ([7]) is replaced by an S'-wave 
bound-state radial function. 
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2.3 Transition probability per time and energy units 



For the Ml transition to the ground state, the gamma width is calculated from 



where is the wave number of the emitted photon. This definition involves 
bound-state wave functions on both sides. 

With the normalization ([8]) of the scattering wave function, the Ml transition 
probability of the process 

^Li(0+) ^ a + ci +7, (10) 



per time and energy units, is given by reduced matrix elements between the 
initial bound state and the final scattering states as (see Appendix A) 

^=^?Sf^^^i(^r.(^)ii-Mfii<)p, (11) 



dE 3h^kad 



where m^r is the nucleon mass. The maximum a + d energy is Q = 2.089 
MeV. The Ml differential gamma width per energy unit to continuum states 
is expressed as 

dr^(0+ -^a + d) _ dW^ 

dE " ' ^ ' 

and the total width is deduced by integration over the energy. 

The Ml operator contains orbital and spin-dependent components. For a gen- 
eral three-body system, it reads, in Jacobi coordinates 
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+ J29s{i)sif,], (13) 
1=1 

where Hn = eh/mNC is the nuclear magneton, Sj are the spins of the three 
particles, and gsii) their gyromagnetic factors. Coefficients A^, Ay and A^y 
are related to the mass and charge numbers as 
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A. 



■X 



Z2AI + ZiAl 
A1A2A12 



A. 



■y 



(Zi + Z2)Al + Z,Al 
AA12A3 
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A^ 



-xy 



I A1A2A3 /Zi 

V AI2A Ui 



Z2 
A2 



) 



(14) 



where ^412 is the reduced mass of the 1 + 2 system; in the present case it is 
denoted as fipn- Variables Px and Py are the momenta associated with the 
Jacobi coordinates x and y, respectively. The matrix elements of the Ml 
operator between hyperspherical functions are given in Appendix B. 



3 Results and discussion 

3. 1 Conditions of the calculations 

3.1.1 Three-body wave functions of the 1+ and 0"^ states 

The initial O"*" wave function is calculated in an a + n + p three-cluster model, 
using hyperspherical coordinates, as explained in Ref. [6]. The same model 
is applied to the ^Li ground state. In both cases, the Coulomb a + p inter- 
action is included, and is taken as a point-sphere potential parameterized as 
ydj) = '2e'^eYi{r / Rc) with a radius Rc = 1-2 fm. Two-body forbidden states 
are removed by using the Orthogonalising Pseudopotential method [7]. The 
central Minnesota interaction ^ describes the n + p system. It is adjusted on 
the deuteron binding energy and reproduces fairly well nucleon-nucleon phase 
shifts at low energies. For the a-\- N nuclear interaction we employ the poten- 
tial of Voronchev et al. [9j, slightly renormalized by a scaling factor (1.008 for 
l"*" and 1.043 for 0"*") to reproduce the experimental energies with respect to 
the three-body threshold (—3.70 MeV for the ground state, and —0.13 MeV 
for the 0"*" state). We truncate the hypermomentum expansion to Kmax = 20 
which ensures a good convergence of the energies. 

The matter r.m.s. radius of the ground state (with 1.4 fm as a radius) is 
found as a/< > = 2.25 fm, a value slightly lower than the experimental 
value (2.32 ± 0.03 fm [10]) (note however that a significantly larger radius, 
2.54 ± 0.03 fm, was found in Ref. [H]). For the excited 0^ level, we find 
V<r^> = 2.56 fm, which is close to the ^He radius. This large value confirms 
the halo structure of this state [12]. The ground state is essentially S = 1 
(96.0%). The ^ = component is 84.4% for ^U{0+) and 82.1% for "^He. 
The ^Li(O^) and ^He hyperradial wave functions are plotted in Fig. [T]for the 
dominant K = 0,2 hypermoments. The overlap (^He|T+/-\/2|^Li(0"'")) = 0.996 
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is in good agreement with the value of Arai et al [12] (0.995). According to 
charge symmetry the short-range parts of the ^He and ^Li(0"'') analog levels 
should be very close to each other. This is confirmed by Fig.[Tl On the contrary, 
the halo components of both wave functions are expected to differ significantly: 
the charges of the halo nucleons are different, and the binding energy of ^Li(0''') 
is much lower. Consequently, the asymptotic decrease of the wave function is 
slower, and matrix elements involving this long-range part should be different 
from their analogs in ^He. 



0.5 




p(fm) 

Fig. 1. Hyperradial wave functions for ^Li(O^) (solid lines) and ^He (dashed lines) 
corresponding to = L = 0, X = 0, 2. The insert shows the same plot in a loga- 
rithmic scale. 

3.1.2 a + d scattering states 

In the foHowing we use four different a + d potentials (the phase shifts can be 
found in Ref. ^): 

(1) The attractive Gaussian potential of Ref. [13] Va contains a forbidden 
state and provides the correct ^Li binding energy [Ei = —1.473 MeV, 
with respect to the a + d threshold). Owing to the presence of a forbidden 
state, it also provides a good fit of the low-energy experimental phase 
shifts. 

(2) To test the influence of the short-range part of the a + d wave functions, 
and in particular of the node location, we also use the potential V^^, 
obtained from a supersymmetric transformation [T3]. The resulting po- 
tential gives the same phase shifts and the same ground-state energy as 
the initial potential, but the forbidden state is removed and the role of 
the Pauli principle is simulated by a short-range core. 

(3) The calculation is complemented by two folding potentials, using the 
deuteron wave function "if a provided by the Minnesota potential. These 
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potentials present one forbidden state. The folding potential V/i is ob- 
tained from the original a + N potential with a renormalization factor 
1.068, which yields the correct binding energy for ^Li; however the quality 
of the S'-wave phase shift is poor. The folding potential Vf2 with a renor- 
malization factor 1.15 describes the S'-wave phase shift accurately, but 
overestimates the binding energy of the ^Li ground state (—2.386 MeV). 

In all cases, the a + d Coulomb potential is chosen as in Ref. [13], i.e. as a 
bare Coulomb potential. 



3.2 Ml properties of bound states 



A test of three-body wave functions is provided by Ml spectroscopic prop- 
erties, which are well known experimentally ^15j . In Table [1], we present the 
calculated values of the magnetic moment and of the B(M1) in ^Li. Separate 
contributions are given for the orbital and spin terms of the Ml operator [see 
Eq. f[T^ ]. In both cases, the contribution of the orbital term is small since 
the dominant component in the ground-state wave function is an S wave. The 
main contribution to the Ml matrix element comes from the spin term. The 
present matrix element corresponds to 5 (Ml) = 7.9 W.u., or = 7.5 eV, 
which are in good agreement with experiment (8.62±0.18 W.u. and 8.19±0.17 
eV, respectively) . The results are also close to those of Kukulin et al. pjj who 
use different variants of a three-body model. 



Table 1 

Contributions (in //^v) of the orbital (L) and spin (S) components to the Ml matrix 
elements. The three-body model is used for the 0^ state. Experimental data are 
taken from |15j . 

(L) (S) Smu Exp. 



Three-body model for l"*" 

/i(^Li) 0.02 0.84 0.86 0.82 

(^-iLill-^f II^SLi) 0-13 2.04 2.17 2.28 
Two-body model for 1+ 



/i(^Li) 0.88 0.88 0.82 

(^-iLill-^f ll^6Li) 0.04 1.53 1.57 2.28 



These matrix elements can also be obtained with a 2-body description of the 
^Li ground state. In that case we use the potential Va to generate the wave 
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functions. Since components with L 7^ are small in the ground-state wave 
function, the two-body model is expected to be a good approximation. The 
l"*" magnetic moment in the two-body model is a simple sum of the proton 
and neutron magnetic moments. In this case, both approaches provide sim- 
ilar results, in good agreement with experiment. However the rms radius in 
the two-body model is a/< > = 2.11 fm, lower than experiment and than 
the three-body value (see Sect. 3. 1.1). In addition, the Ml transition matrix 
element (see Tabled]) provides = 3.9 eV, i.e. an underestimate of the exper- 
imental value. These results suggest that the short-range part of the two-body 
description is too simple. However, transitions to the continuum are more 
sensitive to the long-range part of the a — d wave functions. 

3.3 Ml-transition to the a + d continuum: effective wave functions and their 
integrals 

Since the a + d relative motion is described by S waves, the first and second 
orbital terms of the Ml transition operator f[T^ do not contribute to the 
reduced matrix elements for transitions to the a + d continuum. The orbital 
and spin terms yield nonzero matrix elements only for the ^x = C■y = L = S= l 
and ^x = ^y = L = S = Q components of the three-body wave function, 
respectively. As it will be shown further, the main contribution comes from the 
spin part of the transition operator. The P-wave hyperspherical components 
give small corrections to the process since the 0^ state is essentially S* = 0. 

In order to analyze the 7-decay process to the continuum, we introduce effec- 
tive wave functions and their integrals, in analogy with the /3-decay study of 
the ^He halo nucleus into the a + d continuum We restrict the presentation 
to the dominant spin part. For the initial 0"*" state, let us define the effective 
wave function with hypermomentum K 

uii\R) = {^f^R I dr<P'^{a) ^^ru,{rl (15) 
and the effective integrals 

R 

lf\R)= J dR!uE{R!)u[f{B!), 


K 

where p and a depend on (r, R) , as given in Eq. ([2]) . The normalization factor 
in Eq. (IT^ arises from the Jacobian between the (ii, r) and (a;, y) coordinates. 
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The reduced matrix elements of the Ml operator (spin part) are then directly 
proportional to Ie{R) (see Appendix B). 

0.2 I 1 




Fig. 2. Integrals I^'iR) [Eq.([l6])] at E = I MeV for the q + d potential Va and 
different K values (labels). 

In the following, we analyze two properties: the convergence of the hypermo- 
mentum expansion, and the sensitivity of the effective integrals with respect 
to the a + d potential. Let us start with the influence of -R'max- In Fig- El we 
show the integrals I^\r) calculated at E = 1 MeV with potential Va, for 
different K-values. The dominant contribution at large R values comes from 
the K = 0, 2, 8 components in the ^Li(0"'") wave function. The components 
K = 4 and K = 10 give smaller and comparable effects to the process. The 
contributions of other components are small and not visible at the scale of the 
figure. Similar results were obtained for the (3 decay of ^He jl]. In Ref. [3] it 
was shown that the K = 4 and K = Q contributions are affected by cancel- 
lation effects, which do not occur for i^' > 8. The situation is therefore very 
close to the ^He(0"'") beta decay into the a + d continuum which confirms 
the halo structure of the ®Li(0"'") state, suggested by its large r.m.s. radius. 

In the second step, we analyze the sensitivity of the effective integrals with 
respect to the potential. In Fig.[3l these integrals are shown at E = 1 MeV. The 
potentials Va and Vf2, which provide similar phase shifts and wave functions, 
give results close to each other. This is due to the similar node positions near 5 
fm of the corresponding scattering wave functions. The folding potential Vfi, 
owing to a poor phase-shift description, yields a scattering wave function with 
an inner node shifted to the right (about 0.7 fm), and therefore provides a 
different integral. 

In Fig. [3l we also show as dotted lines, for each potential, the effective in- 
tegrals obtained for the ^He /? decay [3] (notice that in Ref. [1], the factor 
{(A — 2)/A)^^* = 0.74 in the effective wave function was missing). In that 
work, we have shown that a node in the a + d continuum wave functions is 
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0.3 




-0.4 I — 1 

Fig. 3. Effective integrals Ie{R) [Eq. p^ ] at E = 1 MeV for different a+d potentials 
(solid lines). Dotted lines represent the equivalent integrals for the ^He /3 decay [3]. 

responsible for a nearly perfect cancellation effect in the /5-decay matrix el- 
ement. This is illustrated in Fig. [31 for the recommended potential Va-, the 
internal contribution to the matrix element is about —0.30 whereas the exter- 
nal term is about +0.35. The final result is therefore much lower than each 
component individually. This phenomenon yields a strong sensitivity of the (3 
spectrum with respect to the a + d potential. 

Coming back to the 7 decay of the ^Li analog level (solid lines in Fig. [3]), the 
internal parts of the matrix elements are very close to their ^He counterparts 
up to about 10 fm. However, as the long-range parts of the wave functions are 
different in both nuclei, the external contribution to the 7-decay matrix ele- 
ment is significantly larger (about +0.42 for Va)- Consequently, a cancellation 
effect still occurs, but is less important. If we disregard potential V/i which 
does not reproduce the a + d phase shifts, and hence the correct location of 
the nodes in the continuum wave functions, all potentials provide the same 
sign for the matrix element. 



3.4 Ml -transitions to the a + d continuum: transition probabilities 

In Table [2] we give the contributions of different K values to the Ml reduced 
matrix element into the a + d continuum. As we noted above, the orbital and 
spin parts of the Ml transition operator yield nonzero matrix elements only 
with the lx = (!-y = L = S = \ (P-wave) and ix = ^y = L = S = (S'-wave) 
components of the three-body wave function, respectively. As expected from 
the previous analysis, the dominant contributions come from the K = 0,2 
and 8 components. Additionally, the contribution of the orbital part of the 
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Ml transition operator is strongly suppressed (2% at most). 
Table 2 

Contribution of different ^Li(O^) hypermomenta to the Ml reduced matrix elements 
for transitions into the a + d continuum (in IO^'^^at) for the orbital (L) and spin 
(S) terms at several energies 





E = 


0.5 MeV 


E = 


1 MeV 


E = 


1.5 MeV 


K 


(L) 


{S) 


(L) 


(S) 


{L) 


{S) 








-56.9 





-54.2 





-46.4 


2 


0.5 


-85.9 


0.6 


-104.9 


0.4 


-115.6 


4 


3.1 


-36.3 


4.5 


-31.0 


5.2 


-21.0 


6 


1.5 


-11.4 


2.1 


-4.9 


2.2 


1.4 


8 


-1.2 


-51.1 


-1.6 


-61.6 


-1.7 


-61.3 


10 


-0.3 


-22.2 


-0.3 


-23.1 


-0.3 


-20.0 


> 10 


0.3 


-21.2 


0.4 


-14.5 


0.4 


-7.3 


Sum 


3.9 


-285.0 


5.7 


-294.2 


6.2 


-270.2 



To analyze the convergence with respect to the upper bound -Rmax [see Eq. 
f|T6l) ]. we display in Fig. H] the differential width dV^/dE for several values 
of -Rmax (potential Va is used). From Fig. HI one can see that -Rmax = 10 fm 
is far from sufficient. Achieving a precise convergence requires larger values 
(~ 25 — 30 fm), as in the beta-decay calculations of the ^He halo nucleus into 
the a + d continuum [3]. This is not surprising as the halo structure of ^Li(0~'') 
is even more pronounced (see Fig. [T]). 

In Fig. [5l we display the differential width dT^/dE for several a + d poten- 
tials. Contributions from three-body components up to -f^max = 20 are taken 
into account with the maximal relative distance -Rmax = 30 fm. The folding 
potential Vfi shows a picture strongly different from the other ones, with even 
a sharp minimum at about E = 0.8 MeV. This potential gives a poor descrip- 
tion of the a + d phase shift (see Ref. [1]) and hence a shifted node position for 
the a + d scattering wave function. This results in a strong cancellation effect 
as explained in the previous section. The folding potential Vf2 and the deep 
potential Va give close results and the supersymmetric potential V^^ slightly 
overestimates them. 

The integrated 7 widths (from -Z? = up to Q = 2.089 MeV) are given in 
Table [3] for the different potentials. The Gaussian potential Va simultaneously 
reproduces both the ^Li ground state binding energy and the S'-wave phase 
shift at low energies. Additionally, the S'-wave scattering wave function of this 
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Fig. 4. Differential width for Ml transitions into the a + d continuum with the a + d 
potential Va for several values of i?max (in fm)- 
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Fig. 5. Differential width for Ml transitions into the a + d continuum for several 
potentials. 

potential has two nodes at short distances (one due to the ground state, and 
one due to the Pauli forbidden state). The nearly phase-equivalent potential 
Vf2, which also has a forbidden bound state (and hence two nodes at short 
distances) gives similar results. 

The influence of the nodes in the scattering wave function can be tested by 
using potential V^^^ The non-physical ground state of Va is removed by using 
a supersymmetric transformation [T3j. The resulting phase-equivalent poten- 
tial V^^^ has exactly the same ^Li ground-state energy and the same S-wave 
phase shift as Va but its scattering wave functions have one node less at small 
distances. The corresponding width of the Ml transition is about two times 
larger (see Table [3]). In the ^He /?-decay process, this potential strongly over- 
estimates the data [1] . Notice that a very different result is obtained with the 
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folding potential V/i, which has two bound states, but does not reproduce the 
a + d phase shifts and the ^He delayed /3 decay. The shape and magnitude of 
the transition width and probability are strongly different from the result for 

Va. 



Table 3 

Integrated 7 widths for different potentials, and branching ratio BR = T^{0 
a + d)/T^{0~^ — > 1+) (we use the theoretical value r^(0+ — > 1+) = 7.5 eV). 



potential 


(meV) 


BR 




Va 


0.90 


1.2 X 10" 


-4 




0.04 


5.3 X 10" 


-6 




1.08 


1.4 X 10' 


-4 




2.27 


3.0 X 10" 


-4 



Considering the I4 and Vf2 potentials, which are consistent with the data 
on ^He j3 decay, we deduce a recommended branching ratio of 1.3 x 10~^ 
by averaging both values. A previous estimate [2] of the branching ratio 
r^(0+ ^ a + d)/T^{0+ 1+) provides 0.8 x 10"^ This value is close to 
our results obtained with potential Va [13], and is also similar to the branch- 
ing ratio observed in the f3 decay of ^He [17]. Such a branching ratio should 
be observable experimentally. 



4 Conclusions 

In the present work, we have studied the Ml transition process from the 
^Li(0"'") halo state into the a + d continuum and into the ^Li(l"'') ground state. 
Our goal was twofold: (i) to determine the energy distribution of the 7 width 
for the decay into the continuum, and to analyze its sensitivity with respect 
to the a + d potential; (ii) to compare this process with the ^He /3-delayed 
decay. This comparison is a good tool to test charge symmetry in exotic nuclei. 
The ^Li(0+, 1+) states are defined in the three-body hyperspherical formalism. 
The experimental magnetic moment of the ground state and 7 width of the 
0~^,T = 1 state are reproduced with a good accuracy. 

We have shown that the spin-dependent term of the Ml transition operator 
gives the essential part of the matrix elements. In order to test the influence 
of the ^Li bound-state wave functions, we have also used the supersymmetric 
transform [2] instead of the Orthogonalising Pseudopotential method for the 
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removal of forbidden states in the three-body wave functions. The results are 
very similar to the present ones, and were therefore not shown. 

In the Ml transition probability, the K = and K = 2 components of the 
three-body wave function provide about 50% of the matrix elements; conse- 
quently, higher hypermomenta play an important role. The same conclusion 
holds in the ^He jS decay into the a + d continuum, where large K values 
cannot be neglected. 

Ml transitions to the continuum provide a good probe of the halo structure in 
the ®Li(0~'") state. The comparison with the ^He f3 decay shows that the inner 
parts of the matrix elements are very close to each other, as expected from 
charge symmetry. However, the halo parts are different, owing to the different 
binding energies, and different charges of the halo nucleons. In ^Li, the binding 
energy is lower, and therefore the asymptotic decrease of the wave function 
is slower. Consequently the halo contribution is larger in the 7-decay matrix 
element, and even represents the dominant part. This leads to the conclu- 
sion that charge-symmetry breaking is rather strong in these processes. The 
nearly perfect cancellation effect between short-range and halo contributions 
observed in ^He /?-decay is less important here, and the sensitivity with re- 
spect to the potential is therefore weaker. Several a + d potentials were tested. 
The sensitivity is still important (about a factor of 2), but lower than in the 
^He /3-delayed decay. 

The present branching ratio of about 1.3 x 10~^ is consistent with the value 
of Ref. where the authors use a simplified model. The present value is 
based on potential Va which reproduces the ^Li binding energy, the a + d 
low-energy phase shifts, and provides fair results for the ^He /3 decay. It is 
therefore expected to have the same quality for the ^Li 7 decay. An experi- 
mental measurement seems to be possible with current facilities, and would 
provide, in combination with the data on ^He /3 decay, an important step in 
a better understanding of the halo structure in isobaric analog states. 
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A Gamma delayed transition probabilities to continuum states 



Let us assume a bound initial state at energy Ei with spin and parity Jj, vTj of 
a nucleus at rest, decaying to a final unbound state at relative energy E, and 
with spin and parity Jf,TTf. In the final state, both nuclei are characterized 
by spins Ji and I2, and by internal wave functions (j)^'^ and (p^^. According to 
Ref. [is}, the transition probability per time unit is given by 



dW^ = '^-^^dpdPdp^ 5iP + p^)5{E + E,- E,), (A.l) 

where we neglect recoil effects. In ( lA.ip . (p, P,p-y) are the relative, total, and 
photon momenta. The transition matrix element Tji is obtained from 



|T,,r = -^^ E \{^7'''-\p)m^'^''^''^)\', (A.2) 

where (z/i, z/2) are the spin orientations in the exit channel, H!^ is the electro- 
magnetic-emission hamiltonian with polarization q, and is the photon wave 
number. The final state is described by an ingoing wave ^y^^-^ with relat ive 
momentum p = {p,Qp), related to the corresponding outgoing wave '^y^^'^^ 
by 



where K is the time-reversal operator. The outgoing wave function is written 
in a partial wave expansion as 



E {hl2iyiiy2M{iImu\JM)^irimLi^p), (A.4) 

JMnilu 

where J)q^[Qp) are Wigner functions. When the relative coordinate r is large, 
the asymptotic behaviour of the partial wave is given by 



^/f -(E) ^ M^i±lll^^^+iexp(za,) (le{kr) - U'^0,{kr)) 

where cr^ are the Coulomb phase shifts, and and Oi are the ingoing and out- 
going Coulomb functions, respectively. Here and in the following, we assume a 
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single-channel problem or, in other words, that the dimension of the collision 
matrix U is unity. 



After integration over P and p^, Eq. (lA.ip is transformed as 



^2 

dW^ = j^-^\Tf,\^dpdn^. (A.6) 

First, we expand H!^ in electric (cr = E) and magnetic (cr = M) multipoles 
|19] . Then we integrate over the orientations Qp and fi^. We have 



2Ji + l ^ 2A + 1 



^f^lK7/(^)ll-^All^''^')P, (A.7) 



X 



where TVI^ are the multipole operators of order A (coefficients are given, 
for instance, in Ref. [I9]). Let use define 



Using (IA.7P in (IA.6P gives 



dW, _ fik T,{E) 

dE 27r^h^2if + V ^ ■ ' 

where fi is the reduced mass. An interesting case concerns transitions to a 
narrow resonance with energy Er and particle width T. In such the 
scattering wave function can be approximated as ^20j 



hh^ > k ER-E-iT/2 ^ ^ 

where '^'bsa the bound-state approximation of the wave function, and v 
the relative velocity. Using this approximation in (1A.9P and integrating over 
E gives 



W^ = Tf^/h, (A.ll) 
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where T^^^ is the 7 width in the bound-state approximation. This result 
corresponds to the usual definition of the transition probability between two 
bound states. 



B Matrix elements of the Ml transition operator in hyperspherical 
coordinates 



Let us write the three-body wave function ^ as 



(B.l) 



where index 7 stands for (i^^yLS). A reduced matrix element of is obtained 
from 



L J S 
L L' 1 I \ J' L' 1 



x( 



/p, (B.2) 



where we use the notation I = \/2i + 1, and where the integral Ip is defined 

as 



(B.3) 



Matrix elements of iy are obtained by swapping and iy. For the crossed 
term in ffT3l). the calculation is more tedious. We have 



(vl>:;^||^ xpy + yx p.W^';]^,) = 6ssi-f^'^''^'^^'^VQiJyi'J'yLL'J' 



( 0, 



X 



1 i 







1 / 
y ^ 



L J S 
/ \ J' L' I 



^'x ^'y L' 
1 1 1 



^ ^ p I a 1 



(B.4) 



where the angular integral reads 



tt/2 



da cos^ a sin^ a(f)j^ " (a) 



d AL 

+ 



AL 



da tan a cot a 



0^/(a). (B.5) 
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In this expression, A£ = 1 + [^'(^' + 1) — + l)]/2. Integration over a is 
performed numerically. For the hyperradius p, the use of Lagrange functions 
makes the integral very simple. 

For the spin part of the Ml operator, we have 



1/2 



\ J S L 


II 


5" Si S2 


\S' J' \ ^ 







pi 



(B.6) 



where we have assumed that the core spin is zero (53 = 0). 



For transitions to the continuum, the previous formula can still be applied, 
but the final-state wave functions are now defined by Eq. ([5]). It is clear that 
with the restriction to the ^-wave final state, the orbital components and 
ly^^ do not contribute to the Ml transition. The matrix element of the crossed 
term is performed over the Jacobi coordinates. Using the S-wave character of 
the scattering state, we have 



K 



xxyp' 



-5/2 



d 



d 



y dy ^ dx 



(B.7) 



where p and a are given in Eq. ([2]). The spin contribution is obtained with the 
same technique, with the help of Eq. ( IB. 61) . Note that the bra and ket have 
been swapped with respect to Eq. ( ITT]) . The ordering is simply restored with 
a factor — 1/ -\/3. 
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